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ABSTRACT 

We have imaged the M2Iab supergiant star a Ori at 11.7 /um and 17.9 /im using 
the Hale 5-m telescope and an array camera with sub-arcsecond angular resolution. 
Our images reveal the circumstellar dust shell extending out to ~ 5" from the star. 
By fitting two multi-parameter models to the data, we find values for the dust shell 
inner radius, the temperature at this radius, the percentage flux produced by the 
circumstellar material at each of our wavelengths, and the radial dependence of the 
dust mass density. We present profiles of the mean surface brightness as a function of 
radius and the power detected from within annular rings as a function of radius. 



Subject headings: a Ori - stars: circumstellar shells - dust 



1. Introduction 

As stars ascend the asymptotic giant branch (AGB), they shed mass, developing a 
circumstellar shell consisting of atomic gas, molecular gas, and dust. One of the most frequently 
studied AGB stars is the oxygen-rich M2 lab star a Ori, which is relatively easily observed due 
to its high luminousity and proximity, d ~ 130 pc (where 1" = 130 AU). Observations of its 
spatially extended circumstellar shell have been made in many wavelength bands. Mapping at 21 
cm has shown that H I emission extends ~ 1' from the star ( Bowers fc Knapp 1987] ). Large beam 



far-IR spectra have detected [C I] 609 /um in the shell (|Huggins et al. 199^) as well as [O I] 63 



^m and [Si II] 35 /im from what is believed to be the region between the chromosphere and the 
extended shell (|Haas fc Glassgold 1993| ). In the 10 fim band, the shell has been observed using 
both interferometric ( [Bester et al. 199^ ) and drift scan ( |Bloemhof et al. 1984 ) techniques, and 



long-slit 7-14 /xm spectra have revealed spatially extended emission from silicate grains (Sloan et 
al. 19931 ). 



Despite the quantity of previous observations, there are still many unanswered questions 
about the morphology and the optical properties of the dust. The drift scans taken by Bloemhof 
et al. (1984) indicate an asymmetry in the dust shell, implying possible smaller scale structure in 
the circumstellar envelope; however, no fully sampled, sub-arcsecond resolution map of the shell 
has been published. In addition, the size, temperature distribution, and dust condensation radius 
of the silicate grains have not been determined. 
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This paper reports the first successful attempt to spatially resolve the dust cloud surrounding 
a Ori by direct imaging in the mid-IR. We also present several physical models fit to our 
observations. 



2. Observations and Data Reduction 



The observations were made using SpectroCam-10 (Hayward et al. 1993 ) at the Cassegrain 



focus of the Hale 5-m telescope.0 The angular size of the 128x128 Si:As BIB detector's pixels 
(0'.'256) weh samples the diffraction hmit of the telescope ((//S FWHM at 10 /im, (y.'9 FWHM at 
18 /xm), making it possible to resolve the dust shell of the star in the mid-IR. Experience has 
shown that the diffraction limit is frequently achieved with this system. The observations were 
made on 1994 November 13 and 15. Data were taken using 11.7 /im and 17.9 //m filters with 
bandwidths of ~ 1 /im and ~0.5 //m, respectively. The 11.7 /im filter was used for observations 
on both nights, while the 17.9 /im filter was used only on the second. Observations of a Tau were 
interleaved with the observations of a Ori on both nights, in order to monitor the point spread 
function over the course of our observations. An active tip-tilt system was employed, minimizing 
seeing degradation due to telescope movement, and all observations were made with standard 
beam-switched chopping and nodding techniques. 

Individual exposures of a Ori and a Tau were fiat-fielded, coadded, and mosaiced to produce 
images of a Tau and a Ori at each observed wavelength. In the mosaicing process, each input 
image was rebinned by a factor of 3, replacing each pixel with 9 pixels (3x3). In this magnification, 
no interpolation was used; each of the subpixels was assigned a fiux equal to the one ninth of 
the original pixel's fiux value. The images were registered to the nearest subpixel (0'.'083) and 
coadded. The final image of a Ori at 11.7/im consists of 11 coadded frames (total on-source 
integration time of 26.3 s), and the final image of a Ori at 17.9/im consists of 9 coadded frames 
(21.7 s on-source integration). The magnitude of a Tau at each wavelength was determined by 
interpolating between published values listed by Gezari et al. (1993) between 2.2 and 22 /xm, and 
fiux-magnitude relations were derived from Cohen et al. (1992). These were used to fiux calibrate 
the images of a Ori at both wavelengths. 



3. Analysis 

The 11.7 /zm and 17.9 /im data are presented in Figures 1 and 2. Panels (a) and (b) of each 
figure show the a Ori and a Tau mosaics. Figures Ic and 2c display the total fluxes within a series 
of annulli of radius h and width w = 1/12" around the stars, with the a Tau profiles scaled to 



^Observations at the Palomar Observatory were made as part of a continuing collaborative agreement between 
the California Institute of Technology, the Jet Propulsion Laboratory, and Cornell University. 
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match the a Ori photospheric emission. The extra emission from the a Ori shell is clearly seen in 
both plots. 

The mid-infrared emission from a Ori arises from two physically distinct components: the 
stellar photosphere with a diameter of 0'.'044 (Dyck et al. 1992), which is much smaller than our 
diffraction limit, and the extended circumstellar material. Our images represent the superposition 
of these two components as convolved with the instrumental point-spread function (PSF). In 
principle we can generate an image of the shell by subtracting the point source component, but in 
practice it is difficult to determine what fraction of the total flux is due to the point source. As 
an example, Figures Id and 2d show a Ori after subtracting PSF's derived by scaling the a Tau 
images to contain 33% of the total (star + shell) a Ori flux at 11.7 /im and 45% at 17.9 fim. 
(These scaling factors were derived from the spike+ modelling in §3.3.) We found that for 6 < 1 
arcsec, such subtractions were very sensitive to the scaling factor, the stability of the seeing 
(which affected the width of the PSF), and the optical aberrations in the telescope (which caused 
a three-lobed effect that varied slightly between the two stars, especially apparent at 11.7 ;um). 
Because the photosphere emission is so bright, small errors in its removal have drastic effects 
on the remaining faint, extended component. However, at 6 > 1 arcsec, where the photospheric 
contribution is relatively small, the subtraction results were relatively insensitive to the choice of 
scaling factor. 

Both the 11.7 and 17.9 ^um images reveal that at 6 > 1 arcsec the shell flux is distributed 
nearly symmetrically about the star. Unfortunately, the 17.9 jim images of both stars contained 
a filter reflection (at 4% of the flux of the star) to the SE of the star. To reduce the effect of 
this reflection, we did not use the SE quadrant in computing our 17.9 /um radial profiles for the 
modelling described below. 

We first attempted to use maximum entropy and maximum likelihood deconvolution 
techniques in our analysis, but they did not produce stable results for the low surface brightness 
shell because of the extremely bright nearby point source. Therefore, we developed a series of 
multi-parameter models for the shell emission and numerically fitted them to the observational 
data. 



3.1. 3-D Models 

The numerical models we developed for a Ori assume the basic geometry shown in Figure 
3. The star is a a 3600 K blackbody point source. The shell is optically thin to the photospheric 
radiation with dust temperature given by T(r) = To{ro/rY^'^ , where r is the distance from the star 
and To is the temperature at the reference radius To = The mass density of the dust is taken 

as p = p{r); the functional form is dependent upon the model geometry. The dust emissivity e(A) 
is the astronomical silicate function from Draine & Lee 1984. The total flux along a line of sight 
at projected radius b equals the sum of the stellar and integrated dust emission: 
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2hc^ 

h{b)aOri = h{b)phot +P~^ 



p{r)e{X)[exp{hc/ XkT{r)} — i] ^dx 



(1) 



where Xi = , * (2) 

'rf — b'^ if r < ri 



and X is the distance along the hne of sight (see Figure 3). The integration begins at Xj = for 



r > Tj and Xj = y — 6^ for r < (i.e. the value of x at the intersection of the line of sight with 
the inner radius of the shell). The shell starts at an inner radius r = rj, but the integration is 
terminated at x = 8", as continuing the integration beyond this point results in no additional flux 
to 1 part in 10^ in the models considered. In (1), p is the fraction of the total flux due to extended 
emission and Ix{b)phot is the photospheric flux distribution, defined 

^x{b)phot - if 6 > ' 

Performing this integration numerically for an assumed form of p{r), we generated a profile 
of the dust shell emission as a function of the mean radius b. The model was then normalized to 
the shell's total fiux (p/atot) added to the point source {I\{b)phot)- To include observational 
effects, we converted this one-dimensional profile into a two-dimensional image, then convolved 
the image with the a Tau point spread function. Lastly, we extracted azimuthally averaged 
one-dimensional profiles /a(&)cm (CM denotes the convolved model) from the model image for 
comparison with similar profiles extracted from from the original mosaics. 

The models were fitted to the data using a downhill simplex minimization method 
(Press et al. 198^). The simplex method attempts to find the best model by minimizing 



= ^ ~ l\(b)^bs^ d ^' -'■'^i^i^ilyj uncertainty in each point was given by cr = constant, but 
b 

subsequently we tested several other possible definitions of the uncertainty a by allowing a to vary 
as a function of radius (cr = cr{b)). For each of these cases, however, we found that the parameters 
returned from the fits were similar for each of our different definitions of a. We therefore conclude 
that the fitted parameters are robust against both the model details and the boundary conditions. 



3.2. Pov^^er-Lavi^ Models 

Initially, the data were fitted using a simple power-law for the dust density (p oc r~"), with an 
inner radius cutoff. Two models were fitted to the 11.7 and 17.9 /xm data separately, and a third 
model linked the two wavelengths by using common values for the temperature, exponent, and 
inner radius, but separate flux percentages, pu,7 and pn.g. The best-fit values of parameters r^, Ti 
(the temperature at the inner radius of the dust shell), n, and p {pu.j and pir.g for the two-color 
model) for these models are listed in Table 1. 
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The observed and modeled 1-D profiles of the star + shell are plotted in Figure 4. The 

two-color power law model fits the data quite well, with a per degree of freedom value of 1.2 
assuming a constant 5% error on all data points. The power-law model has a serious weakness, 
however: it predicts a value of Cf!Q7 for the inner radius of the shell, in contradiction to values of 
~ 1" indicated by interferometric observations (Danchi et al. 1994). 

3.3. The Spike Models 

To resolve the incompatability of our power law models with previous observations, we tested 
two thin-shell models. The first model, the spike model, assumed a thin shell with thickness of 
only C/.'l. This model did not match the observational data nearly as well as the initial power-law 
model. The value was much larger than that found from the thick shell model (by a factor 
of roughly 100), and the resulting parameter values (sec Table 1), such as rj = 4'.'5, are not 
reasonable, as well as incompatible with the model suggested by Danchi et al (1994). Therefore, 
we rule out the basic spike model. 

The second thin-shell model (hereafter termed the spike+ model) was a combination of 
the spike and extended shell power-law models. We assumed that both the thin and extended 
shells began at the same inner radius, r^, and that the temperatures of the two components 
were equal at this radius. One additional parameter fp was introduced to specify the ratio of 
the density of the thin shell to the density of the extended material {fp = /Othin/pthick), giving 
us p{r) = pthickifp^ifi) + r~'^). This model, like the power-law model, was fitted to the 11.7/im 
and 17.9;um data separately, then was run with the two wavelengths linked by common values for 
To, ri, n, and fp. Numerical results from these models are listed in Table 1 and the azimuthally 
averaged 1-D profiles are included in Figure 4. The spike-|- model also produces a good fit to the 
observed data, with a value of 1.0 per degree of freedom for a uniform 5% error on data points. 
In addition, the model returned inner radii ~ 1", in good agreement with the interferometric 
observations. 

3.4. 2-D Models 

One of the major assumptions made in the model calculations was that the circumstellar dust 
exists in a spherical shell which starts at a well-defined inner radius and extends to a significant 
distance from the star. In order to test this assumption, we also tried a model with a disk 
geometry. Using IRAF software, we measured the eccentricities of both a Ori and a Tau to be 
< 0.05 at 11.7 fim and < 0.1 at 17.9 fim, so a disk would have to be nearly face-on. We again fit 
the disk model at each observed wavelength, and then on both wavelengths simultaneously. Both 
the single- and two-color models were marginally worse than our spherical shell models (Table 2), 
with a value of 1.4 per degree of freedom. This model is included solely for completeness; since 
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a number of evolved stars display spherically symmetric dust shells, and disk distributions are not 
observed in these objects, we believe that the disk model is not a realistic representation. 

4. Results and Conclusions 

We found that the two-color power-law and and the two-color spike-|- models fit the a Ori 
images about equally well (Figure 4). However, because of the disagreement between the power-law 
models and previous interferometric observations, we feel that the spike-|- model better represents 
the a Ori shell. To examine the shape and character of the modeled dust distributions in more 
detail, we created PSF-subtracted images using the flux ratios determined by the models. Plots of 
the mean surface brightness of the dust shell component as a function of radius for both models 
(the two differ in scaling of the subtracted image of a Tau) at both 11.7 and 17.9 /um are shown 
in Figure 5. 

The ratio of the 11.7^m flux to the 17.9/um flux, multiplied by the ratio of the emissivities 
of the silicate dust grains at these two wavelengths (Draine and Lee 1984), indicates the grain 
color temperature as a function of radius. The temperature results for both the power-law and 
spike-|- models are presented in Figure 6. Also shown are the effective color temperature profiles 
of the dust distributions of our best-fit models, calculated numerically using the p(r) density 
distributions and the temperature profile for dust grains in thermal equilibrium, T = To{ro/r)^^^- 

From our two-color spike-|- model fits, we find that the dust shell provides 33% =b 3% of the 
total flux at 11.7 /xm and 45% ± 4% at 17.9 /xm. Prom the value fp, the ratio of the density of the 
thin shell at to the density of the extended material at rj is pthin ^ (75 =t '2^)pthick (we know 
that ptYiin^r — constant, where Ar is the thickness of the thin shell, and are therefore limited by 
the size of our resolution element). Because of the difference in the extent of the two portions of 
the model dust distribution, the thin shell component provides only ~40% of the dust shell flux 
(13% of the total flux) at 11.7/im and only ~30% of the flux (14% of the total flux) at 17.9/Ltm. 

We also find a value for the inner shell radius of = I'.'O it (X'l. This, unlike our simpler 
power-law model, agrees well with previous interferometric observations (Bester et al. 1991). 
While this agreement leads us to believe that this model is a better depiction of the actual dust 
distribution, the dust grain color temperature at rj is only 460 K, well below the condensation 
temperature of ~ 1500 K for silicates (Salpeter 1974). One possible explanation for the low 
temperature is that a stellar wind has pushed the dust inside one arcsecond into the thin shell now 
observed. The dust presumably condensed at a much smaller radius, then cooled as it was driven 
away from the star. This would presuppose that the star has not recently been undergoing heavy 
mass-loss, so that no new dust has condensed near the star. If this scenario is accurate, then we 
can estimate how long it has been since the last major mass-loss event. The last major variation 
in the brightness of a Ori occurred during the period 1941-1945. If this brightness variation 
corresponds to the last major mass-loss event, the mean velocity required to bring the dust from 
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this event to 1" is 12 km s . This velocity is in agreement with the 11 km s~ CO hnewidths. It 
should be mentioned that previous work (Danchi et al. 1994) modeled the circumstellar material 
as a single thin shell (or possible two thin shells), wheras our work indicates that a significant 
portion of the dust emission comes from a very large extended shell. This is not a surprising 
difference, given that the interferometric observations were taken with long baselines, making the 
observations very insensitive to flux coming from large radii. 

While we do not favor the power-law model, it is interesting to note that the power-law 
model color temperature at the inner shell radius (rj = C.'l) corresponds well to the condensation 
temperature for silicate grains. By balancing the flux absorbed by an individual grain with the 
flux emitted by the grain, we find that reasonable grain radii (a = 0.01/im and a = O.l^m) 
produce temperatures {Tdust{fi) = 1740 K and Tciust{fi) = 1890 K respectively) which agree with 
both silicate condensation and our power-law model. 

The temperature of the dust at one arcsec from the star (the inner shell radius) is found to 
be To = 460 it 20K. Examining the spike-|- model plots on Figure 6, we see that this temperature 
is in good agreement with the predicted color temperatures from the p{r) distribution, especially 
for projected radii larger than l". We find an average temperature for the spike-|- model of 205 =b 
35 K, compared to Tmean ~ 300 of Bester et al. (1991). 

The last parameter of the dust shell model which was returned from our minimization method 
was the exponent of the density dependence upon radius for the extended emission, n. The value 
of this exponent from our best-fit spike+ model is n = 0.9 ±0.1. However, n = 2 for a shell 
expanding at constant velocity, and n = 3/2 for a shell expanding at the escape velocity. Therefore, 
any plausible, spherically symmetric, steady state density distribution should have 3/2 < n < 2. 
From this discrepancy, we conclude that the circumstellar dust flow must be episodic. This is not 
surprising, given that the inner radius of the shell also suggests episodic emission. 



The authors would like to thank the staff at Palomar Observatory for their assistance 
and the referee for helpful comments. Figures 1 and 2 can be found in gif format at 
http: / / www.people.cornell.edu / pages / sr22 /Aor: . This research was partly supported by NASA 
Contract 960803. 
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Fig. 1. — Images at 11.7//m of a Ori (a), a Tau (b), and the subtracted image (d) produced with 
the PSF scaled to 39% of the total flux in the a Ori image. Also presented is a plot of the total 
flux observed within annuli of inner radius h (c). Each of the images are 13''3 on a side; the small 
bar in the upper right of each image is 1" long. 



Fig. 2. — Images at 17.9/v,m of a Ori (a), ol Tau (b), and the subtracted image (d) produced with 
the PSF scaled to 49% of the total flux in the a Ori image. Also presented is a plot of the total 
flux observed within annuli of inner radius h (c). Each of the images are 13''3 on a side; the small 
bar in the upper right of each image is 1" long. 



Fig. 3. — A diagram of the physical geometry used for our models. 



Fig. 4. — Plots of the power emitted by annular rings of a Ori (including the star and the shell), 
with the best-fit model overlaid (a) for the power-law model at 11.7 //m, (b) for the power-law 
model at 17.9 /nm, (c) for the spike-|- model at 11 .7 //m and (d) for the spike-|- model at 17.9 /xm 
. Units used in these plots are arbitrary. 



Fig. 5. — Plots of the mean flux from the dust shell, produced by subtracting a scaled PSF from 
the raw image of a. Ori, as a function of distance from the star (a) for 11.7 /7,m and (b) for 17.9 
//m. Also provided is a plot of 25% of the flux of a Tau. Examination of these plots shows that 
at small radii (6 < 1"), the implied shell brightness is highly model-dependent, while at large radii 
(6 > 1") the shell brightness is very similar for both models. 



Fig. 6. — The color temperature of the dust shell, using a two-point blackbody flt from the PSF- 
subtracted images. Also plotted are the predicted effective color temperature curves for each model. 
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Parameter Model 



Spike+ 



Power-law'' 



Disk« 



Spike+ 



11.7 
17.9 

Two- Color 



445±40 K 

510 (-h70/-50) K 

460±20 K 



930±80 K 750 K 

1800 (+1300/-700) K 790 K 
1700±700 K 580 K 



490 K 



11.7 

17.9 

Two- Color 



l^liC'l 
0'.'9±0'.'1 
I'.'OiO'.'l 



0'.'66 ± O'.'l 0'.'58 
0'.'15(-^0'.'2/ - O'.'l) 0'.'42 
0'.'07(+0'.'3/ - O'.'Oe) 0'.'49 



4'.'9 



Pll.7 



P17.9 



f 8 



11.7 
17.9 

Two- Color 
11.7 

Two- Color 
17.9 

Two- Color 

11.7 
17.9 

Two- Color 



0.93 (+0.09/-0.06) 

0.93±0.13 

0.9±0.1 

0.33 (+0.03/-0.02) 
0.33±0.03 

0.41±0.04 
0.45±0.04 

0.986±0.003 

0.976 (-h0.01/-0.06) 

0.986±0.04 



2.0 (+0.2/-0.1) 0.8 

1.4±0.1 0.22 

l.OiO.l 0.01 2.1 

0.35 (+0.03/-0.02) 0.33 

0.39 (-h0.03/-0.05) 0.30 0.05 

0.47 (+0.03/-0.04) 0.40 

0.49±0.05 0.41 0.16 



Two-Color 1.0 



1.2 



1.4 



"These models are included solely for completeness. 
''The temperature of the dust at the inner edge of the shell. 
■^The inner radius of the shell. 
''The density power-law exponent. 

■^The fraction of the total a Ori flux from the dust shell at 11.7/jm. 
^The fraction of the total a Ori flux from the dust shell at 17.9^m. 
^The ratio of the density of the thin shell to the density of the extended shell. 
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